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I. INTRODUCTION 

Decoherence of a quantum system of interest through inter- 
action with an uncontrolled environment is a key obstacle to 
realizing practical quantum information processors. A num- 
ber of methods have been proposed to help deal with this prob- 
lem, including error avoidance methods such as decoherence 
free subspaces (2^4), closed-loop suppression methods such 
as quantum error correction (QEC) [5-7 1, and open-loop sup- 
pression methods such as dynamical decoupling (DD) ifS UTOl 
and the quantum Zeno effect (QZE) |[TTl - [T4l . The typical set- 
ting for the QZE is a sequence of frequent projective measure- 
ments of an observable. When the frequency of measurements 
is high enough, this has the effect of forcing the evolution to 
remain within the eigenspaces of the observable. With appro- 
priate choices of the observable, the QZE can be exploited to 
decouple the system from the environment lfT5l[T6l . 

In spite of the fact that one method uses fast unitary op- 
erations while the other uses frequent measurements, there is 
a conceptual similarity between DD and QZE protection of 
quantum states, in that both are feedback-free methods. In- 
deed, in the "bang-bang" limit of arbitrarily strong and fast 
pulses or measurements, it has been shown that DD and the 
QZE are formally equivalent l!T3l [TTl . However, strong pro- 
jective measurements are an idealization of more realistic, 
generalized measurements, and likewise, real dynamical de- 
coupling pulses are subject to constraints of finite bandwidth. 
It has been shown that DD can work to suppress decoherence 
while allowing for universal quantum computation even with 
pulses constrained by finite width and repetition rate Ifl8lj20l . 
Until recently, an analogous result was lacking for the QZE. 

In our earlier paper [1| we showed how the assumption of 
projective measurements can be relaxed and replaced with 
weak, non-projective, measurements. Weak measurements 
extract less information from the system than the correspond- 
ing projective measurements, and consequently do not fully 
collapse the state 0711221 . In the present work, significantly 
expanding on our earlier paper (TJ, we analyze protocols for 



realizing a QZE using frequent, weak measurements that are 
also non-selective, meaning that outcomes are not recorded, 
with the effect that the state after measurement is the ensem- 
ble average over all of the possible outcomes. We will show 
that these protocols can be used to protect arbitrary, unknown 
states encoded within a stabilizer QEC code, arbitrarily well. 
This will be referred to as the weak measurement quantum 
Zeno effect (WMQZE). Since our protocols involve measur- 
ing the stabilizers of QEC codes — a capability that is taken for 
granted in QEC theory 11231 l24l — but we do not assume that 
we can observe or use the measured syndrome, our assump- 
tions are weaker than those of QEC, and hence the ability to 
perform QEC implies the ability to perform our protocols. 



Weak measurements are in some sense the analog of finite 
bandwidth DD and are both more realistic and more general 
than strong, projective measurements. They capture a large 
variety of experimental imperfections and uncertainties ll25l . 
However, here, as in our earlier analysis, the measurements 
are treated as instantaneous, and a generalization to measure- 
ments of finite duration is still lacking. 



The paper is organized as follows. Section [TT] gives essen- 
tial background on weak measurements, reviews the WMQZE 
protocols introduced in Ref. [ 1 1, and states the main result we 
prove in the paper: Theorem [TJ a distance bound quantifying 
the performance of the protocols. Section [III] describes some 
algebraic structures associated with stabilizer quantum error 
correction codes and the behavior of weak measurements of 
the stabilizer elements with respect to these structures. Sec- 
tion [IV] is concerned with proving Theorem [JJ In Section |V] 
some trade-offs are considered, between the number of mea- 
surements M and the final time r, as well as between M 
and the measurement strength e. Conclusions are presented 



in Section VI Several appendices offer additional supporting 
mathematical details. 
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II. BACKGROUND AND STATEMENT OF MAIN RESULT 



A. Weak Measurements 



Consider a system with Hilbert space Us coupled to a bath 
with Hilbert space Hb- The Hilbert space of the system- 
bath composite is denoted H = Hs ^Hb- Let B(H) de- 
note the space of bounded linear operators on W. A measure- 
ment of the system can generally be expressed as a positive 
operator-valued measure (POVM), comprising a set of "mea- 
surement operators" {Mj£ B(Hs) ®lC acting on 
the system and satisfying the sum rule ^ . M^Mj = 1, 
which map a state g to gj = MjgMj/pj with probability 

Pj = Tr(MjgMj). The observables that are to be mea- 
sured in the encoded WMQZE protocols will be elements of 
the stabilizer group of a given quantum error correcting code 
(QECC). As such, they are unitary involutions, i.e., unitary 
operators S € \J(U S ) <S> 1 C \J(H) such that S 2 = 1, and 
therefore have only two possible outcomes (eigenvalues): ±1. 
The weak measurement of such an observable on a state g may 
be parametrized by the measurement strength e as [25| 

VsM = Ps{e)gP s (e) + P s (-e)gP s (-e) (1) 

where 

P s (e) := a+ (e)P s + a_(e)P_ s , (2a) 

a±(e) = y/(l±tanh(e))/2; (2b) 

and 

P ±s :=i(l±S) (3) 

are orthogonal projections onto the ±1 eigenspaces of S. 
Since -Ps(e) 2 + Ps(— e) 2 = 1> the sum rule is satisfied, and 
Vs,e is a well-defined POVM for every e £ (0, oo). Moreover, 
Vs,e describes a parametrized curve through POVM space that 
interpolates between projective measurement of S at e = oo, 
and no measurement at e = 0. These weak, non-selective 
measurements will form the building blocks of the WMQZE 
protocols to be described herein. It may also be observed that 
this 2- term POVM is unitarily equivalent E3l to the 3 -term 
POVM with measurement operators 



Mi = m 2 = m 3 = VTl 

(4) 

where £ := 2aq_(e)a_(e) = sech(e). This three term POVM 
may be interpreted as a measurement with a particularly sim- 
ple classical error, in which, with probability £, no measure- 
ment takes place, and with probability 1 — ( a projective mea- 
surement of S is performed. It is a completely equivalent 
description of the non-selective weak measurement Vs,e< al- 
though the selective measurements corresponding to these two 
POVMs are not equivalent. 



B. The WMQZE Protocols 

Previous WMQZE work applied mostly to particular states 
E6W291 . with some exceptions ll30l . In order to protect an 
arbitrary, unknown fc-qubit state, as well as to facilitate the 
analysis that is to come and to allow this WMQZE method 
to dovetail easily with other protection schemes like QEC, 
we encode the state into an [[n, k, d]] stabilizer quantum er- 
ror correcting code (QECC) Q [23), with stabilizer group 
S = {Si}f_ , and where So = 1. We assume that the code 
distance d > 2, i.e., the code is at least error-detecting, with 

minimal generating set S = {Si}f =1 C S, where Q = n — k. 
Then every stabilizer element can be uniquely decomposed as 

Si = Yiu=i ^v*"' where rj„ £ {0, 1}, i.e., the stabilizer el- 
ements are given by all possible products of the generators, 
whence Q + 1 = 2^. The encoded initial state go com- 
mutes with all stabilizer elements, and so is supported on the 
simultaneous +1 eigenstate of all the elements of S. For a 
given measurement strength, a weak measurement operator 
Vs.t may be generated for each S G S as in Eq. ([TJ. Since the 
stabilizer group S is abelian, the measurements can be per- 
formed simultaneously, and simultaneous measurement of the 
full stabilizer group can be described by the POVM 

V e := n V S ,e- (5) 
ses 

Similarly, we could measure just the generators S, whence we 
define 

:= II P S ,e- (6) 

ses 

Thus we can now define a weak stabilizer group measure- 
ment protocol in which M evenly-spaced measurements of 
the full group are performed over a time interval [0, r]. The 
state of the system-environment composite then evolves as 

(V e U(r/M)) M (q sb ), 

(V e U(r/M)) M := 
VM (r, T M -i) VM fo-i, Tj- 2 ) ■ ■ ■ P e U{ri, 0), (7) 

where gsB is the initial state of the system-bath composite, 
{ T j = j T /^}fLi are the instants at which the measurements 
are applied, U(t,t') is the unitary evolution superoperator of 
free evolution over [t, if], i.e., 

U(t,t%) = Texp ( [C(s)ds)(-) = U(t,t')(-)tf(t,t'), 

(8) 

where U(t,t') is the solution of the differential equation 
^U(t,t') = -iH(t)U(t,f), with the boundary condition 
U(t,t) = 1, and T denotes time-ordering. Here H{t) is 
the Hamiltonian of the system-bath composite, i.e., H(t) € 
B(Hsb), and the superoperator generator £(t) = —i[H(t), ■]. 
Note that V t in Eq. ((TJ, and more generally S, has non-trivial 
action on the system only. 
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In QEC one measures not the full stabilizer group, but 
rather its generators, in order to extract an error syndrome 
0. It has been recognized that these syndrome measurements 
implement a QZE I31II32I . A weak stabilizer generator mea- 
surement protocol comprises M evenly-spaced measurements 
of the generating set over a time interval [0, t], so that the state 

evolves as (P e U[r/M)) M {g Q ), 

(VMr/M)) M := 
VM (r, t m _ x )V,U fa-i.T^a) ■ • • PeU(n,0). (9) 

This can obviously be an important saving over a full sta- 
bilizer group measurement (P e ). If the measurement is per- 
formed, e.g., by attaching an ancilla for each measured Pauli 
observable (as in a typical fault-tolerant QEC implementation 
11231 ). then this translates into an exponential saving in the 
number of such ancillas. We shall consider both protocols 
in our general development below. 



QS up to a deviation that converges to in the large-M limit: 



C. Distance Bounds 



To quantify the behavior of these protocols, we use a dis- 
tance metric as the figure of merit: 



D[Q S {TlQ%{T)} = -\\e S {T)~Q° s {T)\\ 



(10) 



where the norm is the trace norm (sum of the singular values), 
qs{t) = Tr b[qsb(t)] is the reduced density matrix of the 
system at time r after application of one of the two WMQZE 

protocols, and j?g(r) = T^b[q%b( t )] * s tne ^ na ^ state °f tne 
system under the idealized circumstance that the system re- 
tains its internal Hamiltonian evolution, but does not interact 
with the environment. We shall also require the Schatten oo 
operator norm || -|| (the maximal singular value). 

The principal results in this paper will be the proof and anal- 
ysis of the following upper bound on D[qs(t), Q° s (t)]. 



Theorem 1. Assume an arbitrary pure state g$ = \ips}{il>s\ 
is encoded into an [[n, k, d]] stabilizer QECC with stabilizer 
group S. Assume further that the (possibly time-dependent) 
Hamiltonian H — H$ + Hb + Hsb> where H$ and Hb 
represent the system and bath Hamiltonians and H$b the 
system-bath interaction, is such that = H$ + Hb com- 
mutes with the stabilizer (i.e., is a linear combination of sta- 
bilizer elements and logical operators) and that the inter- 
action term H$b — Xi 9 /i ' s a linear combination of 
error terms detected by the code. Let Jq = \\Hq\\/2 and 
J i = ||-ffss||/2, where Jq and J\ are assumed to be finite. 
Finally, let Q = 2 n ~ k — 1 and q = (Q + l)/2. Then the sta- 
bilizer group measurement protocol {VJJ(t /M)) M protects 



D[Qs(T), Q ° s (T)}<[l + r t (M)] M - 

r_(M) 



[1 + r+(M)] M + T g {M)A + {M)^-\M)- 



r t (M) 

T g {M)A_{M)^-\M) - 



„rJ a 



B, 



(11) 



where the bound B can be expanded in powers ofl/M as 

— — (1 + rJ m ) 



B 



C 9 



1 

M 



O 



1 



, M 2 
where 

c 

f3(M) 

ri(M) 
r 3 (M) 

7± (M) 
A±(M) :-- 



sech(e) 

fri(M) j >Ji 
\r s (M) j <Ji 

1 tJ ( tQJ\ 

+ 1 v e 







Qe 



1 



(12) 

(13a) 
(13b) 

(13c) 
(13d) 



+ 1 

= ~(l + j8 + (l + Q0)C*)± 

^\/(l + /3-(l + 0/3)C 9 ) 2 +4g/3 2 C l (Be) 

QK q (l±+l3) + (l + P)[(l + P)-l T ] 
/?(7± - 7t) 



(13f) 
(13g) 

(13h) 
(13i) 

For the generator measurement protocol (VJJ (r/M)) , re- 
place q by 1 in Eqs. (12 I, (13e), and \13i) . In the strong- 
measurement limit ( e — > oo), both protocols yield the distance 
bound 



J n 

r+ := 
r_ := 



max{ Jo, Ji} 

fri(M) j >Ji 
\r s (M) j q <j, 

'r fl (M) j >Ji 

MM) J <Ji 



D[g s (T),gUr)]<e J ^ 



' _iXL J i T Q \ M 

Qe m + e m \ 
+ 1 J ~ 



(14) 



in. STABILIZER QECCS AND INDUCED STRUCTURES 

A large class of QECCs can be described by the stabilizer 
formalism 871 I231 . which we briefly review. A stabilizer S is 
an Abelian subgroup of the Pauli group G n on n qubits that 
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does not contain the element —I, The Pauli group consists of 
all possible n-fold tensor products of the Pauli matrices a x = 
X, a y — Y a z = Z together with the multiplicative factors 
±1, ±i. All elements of G„ are unitary and either Hermitian 
or skew-Hermitian. Since —I ^ S, all elements of S must be 
unitary and Hermitian, and therefore are involutions [S 2 = 1 
for all S £ S]. The stabilizer code C corresponding to S is 
the subspace of all states \ip) which are left invariant under 
the action of every operator in S (S , |'0) = \4>), VS 1 £ S). The 
stabilizer of a code encoding k logical qubits into n physical 
qubits has Q = n—k generators, and S has Q — 2® elements. 
A set of errors {Ei} in G n is correctable (detectable) by the 
code if and only if EjEj (Ej) anticommutes for all i and j (for 
all j) with at least one element of S, or otherwise belongs to 
S. The normalizer N(S) := {N \ [N, S] = VS G S} C G„ 
is the set of logical operations on the code. 

We now fix a minimal set S = {Si, . . . , Sq} of generators 
of the stabilizer group. This generating set defines a group 
isomorphism B : -> S by B(b 1: . . . , b Q ) = l[f =1 B] 3 , 
where Z2 = {0, 1} is the additive group of integers mod 2. 
The inverse function B^ 1 (S) identifies the subset of genera- 
tors whose product comprises S, namely generator Sj partici- 
pates in the product iff bj — 1. Define for each g £ S a group 
homomorphism a g : S — > Z2 by 

a a {S)~{B- 1 {S),B-\g)) (mod 2), (15) 

where (■, ■) denotes the dot product of the two binary vectors. 
This homomorphism a g then counts (mod 2) the number of 
generators shared by S and g. It is symmetric in that 

a s (g) = a g (S). (16) 

We recall some basic properties of homomorphisms of fi- 
nite groups [33 1 . First, a group homomorphism <f> : G — s- H is 
a map satisfying the property </>(Siff2) = ^(91)^(92) ■ This 
implies, in particular, that </>(1g) = Iff- Both the ker- 
nel if = (jr x {t H ) C G and the image 0(G) C H of 
a homomorphism are subgroups of their respective groups. 
If 0(51) = <j>(g 2 ) then ^(gig^ 1 ) = Iff, so if is 
nonempty for some h G H, then (f)~~ l (h) — Kg for any 
g £ _1 (/i). In other words, all non-empty fibers of 4> are 
cosets of the kernel K. Therefore all non-empty fibers have 
the same cardinality as the kernel, so either |</> -1 (/i)| = \K\ 
or I^O)! = 0. 

Lemma 1. 

a s (g)=0VgeSiffS = t, (17) 

and for S 7^ 1, 

WS 1 m = k/WI = [S|/2 = = g. (18) 

Consequently, as = &s' If an d only if S = S'. So 
{(— 1) CTS (') : S £ S} is the set of all 2® complex irreducible 
representations ofS [34]. 

Proof Since 5(0, . . . ,0) = 1, S" 1 ^) = (0, . . . ,0), which 
yields zero when dotted with any binary vector B^ 1 (g). So 



at{g) = for all g. If S 7^ 1, then S — Sj 1 ■ ■ ■ Sj k for some 
subset of distinct generators S^, . . . , Sj k . Then <7s(Sj\) 
1 for all 1 < i < k. So a s {g) = for all g £ S if and 
only if S = 1. In the case S ^ 1, the image of as is all 
of Z2. Then, the fibers K = erg^O) and cr^ 1 (l) are both 
nonempty and therefore are cosets of the kernel K 1331 and 
partition the group S. Thus, \K\ = Icr^T 1 (0) | = | cr~ 1 ( 1) | and 
|S| = \as\0) U ^(1)1 = 2\K\. Finally, if a s (g) = a s .{g) 
for all g £ S, then 

ass'(g) = Tg(SS') = a g (S)+a g (S') = as(g)+a s ,(g) = 

(19) 

for all g £ S, which, by the arguments above, holds if and 
only if SS' = 1, i.e., if and only if S = S'. □ 



A. Isotypical Decompositions 

The fact that the homomorphisms { ( — 1 ) "a ( • ) : g £ S} are 
the irreducible representations of S ~ leads to a natural 
and well known orthogonal decomposition of the state space 
H into code subspaces Q. 

Lemma 2. The isomorphism B : 1® — > S is a faithful unitary 

representation of in terms of operators on H = Hs ® 
TIb, where TLs — C 2 . There is then a unique isotypical 
decomposition [34] ofH into subspaces 

^ = 0^9' V B = V® a >, a g = 2 k dim(H B ), (20) 

where each V g is an invariant subspace of the representation 
B and the projection of B onto any one-dimensional subspace 
V g thereof is the irreducible representation (— l) CT f( '', i.e., for 
any \ip) £ V g , S\ip) = (-l) " 9 ^^ IV')- Since each V g is one- 
dimensional (because S ~ is abelian), the exponent a g 
is the dimension of the subspace V g . With Tr(l) = 2™ and 
all other elements of S traceless, a g = 2 n ~® dim("HB) = 
2 k dim(Hs) f or oil 9- The 2 k dim(HB)-dimensional sub- 
spaces V g are all orthogonal, and Vt is the subspace stabi- 
lized by S. 

Proof. The isotypical decomposition is a standard result in 
representation theory, following from Schur's lemma [34). If 
|V>) G V g , then S\ip) = (-l) a ^\i;) because (-l) a ^ is the 
irreducible representation associated to V g . By ll34l Corollary 
2.16], 

a 9 = 4 E Tr(S)(-ir* (S) = 2-<3 dim(H B ) 
ses 

= 2 k dim(H B ) (21) 

using the fact that S = 1 has Tr(l) = 2", and all other 
S £ S have zero trace. If \ip g ) £ V g and |t/)^) G Vh, then 
since S c U(W), 

= (-l)°"« h(s) (V fl |^h) (22) 
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for all S £ S. If gh 7^ 1, i.e., if g ^ h, then by Lemma [T] defined by 
\ a gh(l)\ = l S l/ 2 ' so mere exists S £ S such that 0^(5) = 
1. Therefore (ijj g \iph) must be zero and the subspaces V g form 
an orthogonal decomposition of H. □ 



(26) 



ses 



In the language of QEC [7], each of the subspaces V g can 
be thought of as encoding k qubits, but only Vt is stabi- 
lized by S [i.e., S\i>) = (-1) CT ^ S )|V> = \ip) V5 e S iff 
a g (S) =0V5eS <=> 3 = 1 by Lemma [l). Hence V t is 
typically chosen as the stabilizer QECC. With this choice, the 
remaining isotypical subspaces are interpreted as "syndrome" 
subspaces, where g labels the syndrome. Namely, after an 
error that is detectable by the code takes place, it maps V\ to 
one of the other subspaces V g . A measurement of all the stabi- 
lizer generators reveals the label g, in that g = B{b\ , . . . , bg), 
where bj G Z2 is (1) if the measurement of generator Sj 
yielded eigenvalue +1 (—1), with j G {1, . . . , Q}. 

Using the group action of conjugation by S, the space B(H) 
of all linear operators (complex matrices) on T~L = Hs ® "Hb 
and the space Ai(T~L) of Hermitian operators may be similarly 
decomposed into isotypical subspaces indexed by S which are 
orthogonal under any inner product invariant under conjuga- 
tion by S (such as the Hilbert- Schmidt inner product). Then 



B(n) = ($W g z M(H) = ($W g , 
g es g es 



(23) 



where, for any g G S, the subspace (respectively W g ) is 
the space of all operators (resp. Hermitian operators) A g with 
the defining property that they satisfy SA g S = (-lY^^ A g 
for all S G S. From this description, it is clear that W g C 
Wg" for all g £ S. For more on these decompositions, see 
Appendix [A] 



Lemma 3. The isotypical decompositions in Eq. \Ti\ impart 
to B(H) the structure of an S-graded associative algebra 4351/ 
and to M{'H) the structure of an S-graded Lie algebra ( under 
the Lie bracket A, B <-> i[A, B]). 

Proof. For any g, h G S and any A g £ W g and A h £ W^, 

SA g A h S = SA g SSA h S = (~l) a ^ S)+a " {S) A g A h 

= (-iy° his) A g A h , (24) 

so that A g A h G Wf h , and therefore W^W^ C Wf h . Like- 
wise for any A g G W g and Ah £ Wh, 

S(i[A g ,A h ])S = i[SA g S, SA h S] 

= (-ly^^+^hiAg^Ah] 

= (-iyMS) i [A g ,A h ], (25) 



so that i[A g ,A h ] G W gh , and therefore i[W g ,W h ) C W, 



gh' 

□ 



Finally, we can define the orthogonal projections into these 
isotypical subspaces as follows. 



is the orthogonal projection into the subspace Wg. Restricted 
to AA(H), this same operator defines the orthogonal projec- 



tion into W g . 



Proof. For any A G B(H), and any g, S' £ S 

s'(v g (A))s' = -j- j2(~ i y AS)s ' sAs ' s 



ses 



±-s2(-ir^ s ' s) sAs 



ses 



ses 



= (-1) CT ' (S > S (A), 
so that V g (A) £ Wg 2 . Moreover, we find 

Y,rM) = ^E fa- 1 )^) 

ges 1 1 ses \ g es J 

since by Lemma [l] E<,es(- 1 )' T9(S) = for S + 1 
equals |S| when S = 1. Therefore, since the subspaces WS 
are mutually orthogonal, the operators V g are orthogonal pro- 
jections. Finally observe that if A is Hermitian, then V g (A) 
is Hermitian as well, so V g also describes the orthogonal pro- 



(27) 



SAS = A, (28) 



and 



jection M(H) >-)• W g 



□ 



Note that Eq. ( |2"7| ) shows that V g (A) coincides with the 
defining property of the operators A g (Hermitian or not), so 
that we can equivalently define A g :— V g (A). 



B. Measurement Operators 

With respect to these isotypical decompositions, the actions 
of the measurement operators defined in Section [i] take par- 
ticularly simple forms. Recalling Eqs. ([T| and ( |2a] i, for any 
S £ S and e > 0, the effect of the weak measurement of the 
stabilizer S is given by 

VsM) ■= Ps(t)AP s (e) + P s (-e)AP s (-e) (29) 
= E E a sAbe)a S2 {bt)P SlS AP S2 s 

b=±l Sl,S2=± 

for any A £ B(H). Then 
Lemma 5. For any g, S £ S, any A g £ W^, and any e > 0, 



Lemma 4. For any g £ S, the operator V g : B(Ti) 



rsAAg) = C si9) Ag, 



(30) 
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where £ = sech(e), and therefore, using Eq. ((28), 

ses se^o) se^^i) 

=^S,oo + C(l - V S ,oo) = (1 - C)^S,co + CI- 



(31) 



Proof. Observe first that, from the definition of W g , 

1 

2 



A g P S2 s = \Mt + S2 s) = \(\ + S2 (-iy^s)A 



g 

(32) 



so that, using the facts that P| s = P± s , PsP-s = P-sPs = 
[in Eq. ((33c)], P s + P_ s = 1, a\ (e) + a| (-e) = 1, and 
Q + (e)a_(e) + e)a_(— e) = scch(e) = £ [in Eq. ( (33e) ], 



b=±l si,s 2 =± 

-PsiS-P S2 („l)"s(s>sA? 



(33a) 



= E E a si( be ) a S 2(-iysM{be)x 

b=±l St,S 2 =± 

P SlS P S2 sA g (33b) 

= E Y as ( be ^ a s(-ir s{a) ( beSlPssA 9 ( 33c ) 

b=±l s=± 

if <r s (g) = 

E&=±1 O) E s=± ^S^ s 

if cts(s) = 1 



A s a s (g) = 
<Ms) = 1 



(33d) 
(33e) 



which proves Eq. ((30). Equation ( fJT) follows from the obser- 
vation that E g Go- _1 (o) ^3 * s * ne projection into the subspace 
of operators that commute with S (the commutant or central- 
izer of 5), which is precisely the strong (von Neumann) non- 
selective measurement of S, Vs,oo- Q 

We are now ready to see the effect of the complete POVM 
defined in Eqs. |5) and |6). 

Lemma 6. For any g g S, A g £ Wg, and e > 0, the weak 
measurement V e of the full stabilizer group has the effect 



A 9 g = t 
C q A g else. 



(34) 



so that the weak measurement may be written as 



v e =Vt + ( q Y. v 9 = A + - A) 
=(i-C q )A + C q t = (i-c q )r 00 + cn. (35) 



Proof. Since P g P h = for g ^ h, 

^ = n -Ps* =iie c s(9s) v gs = e n ^ {s) a 

Ses sesgses gesses 

= ^ C k 9 - 1 Wlp 9 , (36 ) 
ses 

and from Lemma [T] 



kr 1 (i)l = 





q else. 



(37) 



which proves Eq. ( (34) . Equation ( (35) follows from the same 
reasoning as used for Eq. ( (3T) in the previous lemma. □ 

Lemma 7. For any j e S, 4 9 6 Wj", ana? e > 0, f/ie weak 
measurement V e of the generators of the stabilizer group has 
the effect 



V e (A a ) = C lBn ^ 1{1)l A 9 



(38) 



so that the generators-only weak measurement may be written 
as 



E 



Va 



c=1 9 e{hes ■. (i)\=c} 
Q 

= ^+EC C E A jr -S ic - (39) 

C=l l<Jl<-<ic<<3 

Proo/ Since P g -p h = for .g ^ h, 

= n ^ = n e ^ ^ e n ^a 

ses sesgses ses ses 

= J2C l§n <^V g , (40) 
ses 

and from Lemma [T] and the definition of o~ g , 



^'Mil o) Li 1 <«> 



□ 

We note that the result of Lemma(7]is stronger than that re- 
ported in our earlier work [ 1 ], where we used the lower bound 
1 for g 7^ 1 in place of Eq. ( f4T) . 



IV. ANALYSIS OF THE DISTANCE UPPER BOUND 

In this section, we analyze the behavior of the distance 
D[qs{ t )t Q%{t)] [Eq. ( (TO) ] and show that it converges to 
in the limit of large numbers of measurements. The Hamil- 
tonian is orthogonally decomposed as H(t) = J2 g es -^s(*)> 
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where H g (t) is the component of H(t) lying in the isotypical 
(error syndrome) subspace W g . This yields the superoperators 
C g = i[H g , ■]. We also write H = H t and Hi = X) s /i -^3' 
and define Jo = 2|ji?o||oo and Ji = 2||i/i|| oc , where || • Hoc 
denotes the norm, i.e., ||iJo||oo = ess sup tg [ 0T ] ||iJo(*)|| 
in which ess sup denotes the essential supremum and || • || de- 
notes the Schatten oo norm, i.e. the maximum singular value. 
This guarantees that ||£i||oo < Jo and ||£ 9 ||oo < Ji for all 
g 7^ 1 in S. These finite bound conditions may also be shown 
to imply rapid decay of the noise spectrum at high frequen- 
cies, guaranteeing an effective spectral cutoff; conversely, an 
insufficiently rapidly decaying noise spectrum implies that our 
finite bound conditions are not satisfied (see AppendixFE}. 



it is found that 

, . „, /Af-1 \ ,,/M-2 M-l 
Qsb(t) = VM [ -^^t,t) VM f -^^t, -^^t 



M 



M ' M 



■ ■ ■ VM |^0, — tJ q sb = Qsb + 

M oo 

E E E E vi^M-^ 1 )- 



sj es 3 
lliS-M 3=1 n 



. . . r p % ^y ri 



X 



-i^iQSB) (46a) 

M oo 

-QSB+J2 E E E c 

77=1 ii,...,i 7) =l ?eNi 5 j es ! j 

||i|ll<Af 3=1,...,, 

4 1 /- il+ i(« 1 )---4vi- il -.^(« r ')(^B), 

(46b) 



Let N = {1,2,3, . ..} and N = {0, 1, 2, ... } and observe 
that the unitary superoperator ((8) describing the joint system- 
bath evolution between successive measurements can be writ- 



ten as 



U 



T, T = 1 

M M 



where fl € {0,l} n , fij = w^^-i • • • g } ), and g = 
a\ ■ ■ ■ ol\ . . Then 

r>fe(r),^(r)] = \ ||Tr B [g SB (r) - 0° B (r)] ^ 



(47) 



£(ii)£(t 2 )di 2 dti 



E E 

fe=0 a£S k 



(42a) 
(42b) 



by the Dyson expansion and by the isotypical decomposition 
C = J2 g es £g> wnere S fc is the set of all fc-tuples of stabilizer 
group elements, £° = 1, and for any k > and any a £ S k , 



4(5) := f£[ ftk T ■ ■ ■ jhl^Ati^aAh) 
■ ■ ■ C ak (t k ) dtfc • • • dt 2 dti. 
Lemmas [3] and [6]imply that 



(43) 



p: i 4 + i- il -...- ij (^)-"4 + i- <1 -...-i,(^)(te B ) = 

(44) 

where ^ = v[(a{ ■ ■ ■ of.) • • • (a? • • • afj] and u : S -> 
{0, 1} is defined by 



v (g) = for g ^ 1 and = 1. 



(45) 



This follows from the fact, implied by Lemma [3] that a com- 
position of Hamiltonian superoperators C ai ■ ■ ■ C ak will map 
the initial density matrix (an element of the isotypical space 
W%) to the isotypical space W ai ... ak , and that this space, by 
Lemma|6j determines the action of the measurement V e . Then 



where 



W := Tr, 



M oo 

E E E E c-(-): 



rj— 1 h,...,l Y) = 1 ie^ 



II 1 S M j — 1,. . . ,7] where ,3^0 



(48) 



and 



S := Tr, 



M oo 

ew+E E E E 

II i 111 < M j=l,...,rj where ^=0 



~U°(r)(g SB )] , 



(49) 



with U°(t) the system-bath unitary evolution superoperator 
generated solely by Hq. 



A. The "Weak" Term 

The "strong" term S represents the strong-measurement 
limit of the norm argument, i.e., the limit as e — > oo. The 
behavior of ||<S|| i will be analyzed in Section IV C and shown 
to vanish as M — > oo. In the remainder of this section, we 
study the behavior of the "weak" term 1 1 W 1 1 i - To that end, 
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\r ii (r/M)fc 



observe that for any 1 < j < M and any a G S , 

\\^^)\\<\\C ai \\oo---\\C ak \\ cc J^ r jt lT --- (50) 
J*-! 1 dtfc • • • di 2 dti = ||£ Ql ||oo 
where ||£ Q Joo := ess supg< t < r ||>C ai (t) |. so that 

M oo 

iiwiii<E E E E c fl 

ll*lll< M j — 1,..., ?7 where^O 



II^m-u+i^ 1 )!!---!!^!^-.^^) 



(51a) 



< 



E E E c 

^ = 1 i £N i7 K{0,1}1 
1 1 Til , <M 



E 



(r/M)' 1+ -+'" 



e n^ (ft)' 

(51b) 



where 71(g) is defined by 

S 36(0,. ...Q)' ll-^Qllloo ' ' ' HAjJU < ll(g) 

S ai S ai =g 

eU ( l s )J$J l r 3 fis(9), 



(52) 



in which fi (g) denotes the number of ways of generating g £ 
S from I non-identity group elements (i.e., g = S ai ■ ■ ■ S ai ). 
We have used in Eq. ( |5T| the triangle inequality and submul- 
tiplicativity 



\AB\\ < \\A\\\\B\\, 



(53) 



(for any pair of operators A and B) and the fact that if g is a 
normalized state then ||g||i = 1. We have also used the norm 
inequalities 

||AB||i < ||i4||||B||i and \\Tr B [A] || x < \\A\\i, (54) 
valid for any operators A and B acting on "Hs ® ll36l[37l . 

To find a closed expression for fi(g), first note that fi(g) is 
constant over all g 7^ 1. Then for any g = S ai • • • S ai 7^ 1 
it may be observed that g 1 = S ai ■ ■ ■ S oll _ 1 can be anything 
other than g (otherwise S ai = 1, which is forbidden) and 
then g' and g uniquely define S ai . One such choice of g' is 
g' = 1, the other Q — 1 choices are non-identity. The same 
logic applied to the case g = 1 shows that all Q possible 
choices for g' are non-identity. It follows that the quantity 
fi(g) satisfies 

fl(9) = E fl-^9') = (Q~ l) f l-i(9) + fi-iW (55a) 
fl(t)=J2^9')=Qfl-l(9)- (55b) 
Therefore /; (g) satisfies the linear recurrence (see Appendix 



fi(g) = (Q-l)fi-i(g) + Qh-2(g). 



(56) 



The characteristic polynomial of this recurrence, x 2 — (Q — 
l)x — Q = 0, has roots —1 and Q, so fi(g) has the general 
form fi(g) = AQ l + B(—l) 1 . It is easily seen that f (g) = 
and fi (g) = 1, so A and B may be found by solving the linear 
system 



which yields 



so that 



Q + l 



A + B = 
AQ-B = 1, 



B 



(57a) 
(57b) 



o l - f-iV 

/*(g)= „ , forget f t (t) 



Q + l 

Q i +Q(-iy 



(58) 



Q + l 



Q + l 



Now, since 



s=0 



E J V{- s a'- s -(J + aJi)' 



(59) 



(60) 



for any fixed a € K, it is readily seen that 

f v (Jo + QJi)' - (Jo - Ji) ; 
Tils') = 7^—, forg^l (61a) 



71 OO 
Defining 



Q + l 

(J + QJi)' + Q(Jo- Ji)' 
Q + l 

~ (r/M)S(ff) 



r 9 (r/Af):=E 



z=i 



/! 



(61b) 



(62) 



we have 



r 9 (r/M) = e- 
r t (r/M) = 



e M -e~TT 

Q+l 



e »f 



TfjJl T.7] 

e M +Qe M 



for 3 ^ 1 (63a) 
1. (63b) 



Observe that for g 7^ 1, 

iVr/M) - r 9 (r/M) = e 11 ^ - 1 

J > for all t > when J > Jx 
I < for all t > when J < J x . 



(64) 



In addition, the cardinality of the set of 77-tuples (gi, ■ ■ ■ , g v ) 
such that v(g v ■ ■■ 9j) = Mj for all j is QH^H 1 , because gj must 
equal <?,,••■ gj+i when /Uj = 0, and gj may be any of the Q 



9 



other elements of S when fj,j = 1. Therefore, 

\ll+-+ln V 



E 



(r/MY 



e n-Wifo) 



;.!... 7 ! 

Ji,...,/„=1 3i,...,s„es j- = i 



< 



e n r *( r / M ) ( 65a > 

' Q\mir g (T/M)[T 1 (T/M)] r >- 1 J > J, 
Q^[T g (T/M)] v J <Ji, 



(65b) 



since each product in the sum contains at least one T g (g ^ 1). 
Let 



M 



0(M) := ^^(r/Af)"- 1 ^ ^ ^qIImIU (66) 



7,=1 

O + i 



(■£{0,11' TeNI 
||?||i<M 



/?(r/M) := 

' v i tJq / tQJ-i 

T % {t/M) = feTT 



Qe 



r 9 (r/M) = ^e^ (e^-e-*) 
Therefore we have 

|| W||i < T g (T/M)cj>(M). 



(67) 
(68) 

- 1 Jo > Ji 
Jo < Ji- 



(69) 



reasoning, the number of ordered fc -partitions of n or less is 
(?) , seen by considering the placing of k separators between 
n + 1 objects, yielding ji + ■ ■ ■ + jk+i = n + 1, and then 
discarding jk+i- So 



E E c^Q 11 * = 



fl^o ||?|| 1 <M 



n , N M-(rj-r) 



Eo r 



E r 



u- 1\ (M-u 



r — 1 J \ rj — r 



(71) 



and we find that <fi(M) in Eq. ( |66] i is equal to 



'/ 



^(Af) = ^/3(r/M)"- 1 ^Q 
M—(ri-r) 

E e 



X 



r=l 

it - 1\ /M - u 



r — 1 J \ rj — r 



(72a) 



,1/ 



r},u,r=X 



r J \r — 1 / \ r\ — r 

(72b) 



where the limits have been extended to 1, . . . , M since the 
additional terms are all zeros. 



A closed form will now be derived for T g (f>(M) and shown 
to converge to zero as M" 1 for arbitrary fixed r > and 
£ £ (0, 1) as M — > oo. To that end, let r represent the number 
of non-zero elements in the vector, and u represent jl ■ i. 
Then 



E E e*«" 



ME{0,1}1 i €m 
fi* \\~i\\l<M 



M— (77— r) 



2_j C u #{ or dered r-partitions of u} : 



^{orderered (77 — repartitions of M — u or less}, (70) 

where an ordered fc -partition of n is an ordered set of fc pos- 
itive integers that sum to n: ji H + jk = n. The number 

of such partitions is as may be seen by considering the 

placing of k — 1 physical separators between a linear arrange- 
ment of n physical objects (see Fig.[T]l. Moreover, by the same 



FIG. 1. Example partition of six elements into three sets by choosing 
the positions of two separators from among the five possible gaps 
between adjacent elements. 



B. Linear Recurrences 



Let <E> be the summand of Eq. ( f72| ), i.e., 

$(M,u,r,,r) -^P^CQ 



7,-l jC u n r/ ? A (U-1\(M -U 



r-1 



r\ — r 



(73) 



where we regard /3 as a fixed quantity for the moment (its M 
dependence will be reintroduced later). It may be observed 
that $ is hypergeometric in all four variables, i.e., $(M + 
1, u, 77, r) /$(M, u, T], r) is a rational function of M, u, 77, and 
r; likewise for the corresponding shift ratios on u, 77, and r. 
As a result, <£> admits a linear recurrence relation such that the 
coefficients are independent of u, 77, and r, and the recurrence 
may be found algorithmically [38 1. Wegschaider's MultiSum 
package [39, 40] for Mathematica, implementing a version of 
Sister Celine's method [38], yielded the following telescoping 
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linear recurrence for this summand: 



with roots 



Z(l+P + Qp)$(M-2,u,r ] ,r)- 

(1 + p + (1 + Q/3)f )*(M - 1, u, 77, r) + $(M, «, 77, r) 



/3$(M - 1, u, 77, r + 1) + $(M - 1, u, 77 + 1, r + 1) 



$(M, u, 77+ l,r+ 1) 



A, 



f$(M- 2,u,?/,r + l) + (1 + £)x 



$(M - 1, u, 77, r + 1) - $(M, m, 77, r + 1) 



(1 + 0)£$(M - 2, «, 77, r) + (1 + P + x 



$(M - 1, w, 77, r) - $(M, u, 77, r) 



(74) 



where A r) is the forward shift difference operator on 
the variable 77: A v X(M,u,r],r) = X(M,u,rj + l,r) - 
X(M, u, rj, r) for any expression X; likewise for A r and A u . 
The recurrence ( f74] i may be verified by expanding the A ex- 
pressions and collecting terms, leading to the equivalent equa- 
tion 

= Qp£$(M - 2, u, 77, r) + p£$(M - 2, u, 77, r + 1) 
+ £$(M - 2, u, 77 + 1, r + 1) - Q/3£$(Af - 1, u, 77, r) 

- £$(M - 1, 7i, 77 + 1, r + 1) 

- /3$(M - 1, tj + 1, 77, r + 1) 

- $(M - 1, u + 1, 77 + 1, r + 1) 

+ $(M,u+l,77+l,r + l). (75) 

Dividing Eq. ( |75] > by Q$(M — 2, u, 77, r) and multiplying by 
the common denominator r(r+l)(M— w — 77+r— 1), Eq. ( |74| i 
is then found to be equivalent to 

= f3£r(r + 1)(M - it - 77 + r - 1) + ££(77 - r) 2 (ii - r) 
+ /3£(?7 + l)(u - r){M -u-rj + r-1) 

- /3£[r(r + 1) + (77 + l)(w - r)](M - u - 1) 

- ,S£(»7 ~ ^) 2 u ~ P€(v + l)u(M -u-rj + r-l) 

+ PZ(rj + l)u(M-u- 1), (76) 

which is easily verified to be true. 

Now we wish to turn this recurrence on the summand 
$(M, u, 77, r) into a recurrence on the sum <fi(M). Do- 
ing this requires summing both sides of Eq. f74"] > and using 
the fact that the right hand side leads to telescoping sums: 
5Zfc =1 AfeX(fc) = X(n + 1) — X(l). This computation, car- 
ried out in Appendix [C] results in the inhomogeneous linear 
recurrence relation (with constant coefficients) 

(1 + P + QP)&{M - 2) - (1 + p + (1 + QP)£)<j>(M - 1) 

(77) 



+ 4>{M) = QP{l + p) M - 2 t 



This recurrence for the sum has the characteristic polynomial 

x 2 - (1 + P + (1 + Qp)£)x + (1 + /3 + Q/3K = (78) 



7± 



(1 + /3 + (1 + QP)i) 



{I + p + (1 + Q0)f) - 4(1 + + 
± s (79a) 

(l + /3+(l + Q/3)C) 



± 



(l + /3-(l + Q/?K) +4Q/?2£ 



(79b) 



Furthermore, by plugging a sequence of the form <p(M) 
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is 



a(l + /3) M into Eq. (|77}, it may be seen that 
a particular solution to the inhomogeneous recurrence. Thus 
Eq. ( |77| ) admits the general solution 



<f>(M)=A+j?- 1 +A^ r 



M-l 



(1 + /3) 



M 



The initial values of 4>(M) may be worked out to be 

0(1) = 

0(2) = Qf(l + 2/3 + (l + Q/3)e). 



(80) 



(81a) 
(81b) 



Thus, we seek A+ and A_ such that A + + A- = Q£ + (1 + 
p)/p and A +7+ + A_ 7 _ = Q£(l + 2p + (1 + Q/3)£) + (1 + 
P) 2 /P- Solving this linear system for A + and A_ yields 

, Q^(7±+/3) + (1 + /9)[(1 + /3)-7t] rs9 , 

A± - ^ r . (82) 

With these values for A + and A_, and now regarding /?, 7 ±, 
and A-t all as functions of M, Eq. ( |80"] l is an exact closed 
expression for the sum cj)(M). So we finally conclude that 



II < 



r ff (M) 



/3(M)A + (A/) 7 f- 1 (M)+ 



/9(M) 

/3(M)A_(Af) 7 5" 1 (M) - (1 + /3(Af)) M 



where 



/3(M) 



MM) j >Ji 
MM) j <Ji. 



(83) 



(84) 



C. The "Strong" Term 



Following arguments similar to those elaborated thus far, 
we may now readily compute an upper bound for the "strong" 
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term ||<S||i and examine its asymptotic behavior. Explicitly, 



M 



\\S\\i 



E E E 

77=1 /l,...,i^=l ?EN»> S36S ! 3 

ll«lll<*f fJ = Obulnolallc,{=I 



(85a) 



< 



E E E E 



77=1 i 1 ,...,i I) =l iefi>7 



Hi-*' £=0tatm>tallc,{=l 



\^ 1+1 &)\\>>-\\z l & +1 _ il _..._^cr')\ 



(85b) 



< 



E E E 

»)=1 ii,...,i„=l *SN»! 

||?||l<M 



(r/M)' 1+ -+'" 
Zi! ■■■/„! 



n^(i)-Jo i+ '" + '' 



EC!') [n(r/W)-(^-l 



[l + ri(r/M)] M -e 



A/ 



g + i 



(85c) 

(85d) 
(85e) 
(85f) 



we obtain, using Eqs. ( |83j > and ( |85e| i, the upper bound on the 
distance metric 

D[q S (t) ) q s {t)]<^(\\S\\ 1 + \\W\\i) 



< 



[l + r^M)] 



m r_(M) 



[i + r + (Af)] M - e 



I\(M) 

+ r fl (M)A+(M)7f -^M) + r 9 (A/) J 4_(Af) 7 M - 1 (M) 

(88a) 

(88b) 



OO 1 



where 



-^-( l + T J 



1 - C 9 



(89) 



This completes the proof of Theorem [T| T he detailed asymp- 
totic analysis is carried out in Appendix[D| It may be observed 
that the rate of convergence (i.e., the coefficient of M^ 1 ) 
blows up as t —> 00 and also as £ — > 1 (i.e., the weak mea- 
surement limit as £ -> 0). 

In the next section we analyze the conditions under which 
the distance bound converges to zero. 



V. PARAMETER TRADE-OFFS 



A. Trade-off between r and M 



D. Total Distance Upper Bound 



Combining the results for the "weak" and "strong" terms, 
and recalling that A± (M) and j± (M) depend on the form of 

/?(M) 



Ti(M) 



P{M) = { 



r s (M)=eT 



rg£i Tjf 

e m +Qe ~fT~ 
Q+l 



tQJ 1 _ t.J 1 

e m — e m 
Q+l 



1 J > Jl 



Jo < Jl, 
(86) 



with q = |S|/2 = (Q + l)/2, J m = max{J , Jl}, and 
fr t (M) J >Ji _ 



r 



T g (M) J <Ji 



Ft(g) Jo > Ji 
T t (M) J <Ji, 
(87) 



The asymptotic behavior of the distance upper bound, pre- 
sented in Section IV D and derived in Appendix [D] shows that 
the bound converges to zero as 1/M when the time r is held 
fixed. Here we consider how the number of measurements 
needed for convergence depends upon the final time r. In or- 
der to address this question, we consider the convergence of 
the bound when r is allowed to increase with M. If r in- 
creases too quickly, the distance upper bound will no longer 
converge, so we seek the fastest growing function f(M) such 
that when r = f(M), the bound still converges to zero. 



Lemma 8. Allowing only r and M to vary, the distance upper 
bound (|88b|> converges to zero as M — > 00 provided 



tJq = alog(M) 



a < 1 ifJo> Ji 
a<% if Jo < Ji ' 



Proof. Suppose that tJq — alog(M) for some a > 0. Then 
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defining A = Ji/Jq, 



r 1 (M) = e 



/ QAalog(AJ) Aa log(M) 

aios(M) / e h + ye M 



= P M 



V o+i 

alog(M) (1 + QA 2 )a 2 log 2 (M) 



M 2M 2 
log 3 (M) 
M 3 



(91a) 



r s (M) = e - 



/ QXa log(A/) Aa log( A/ ) 

g iog(M) / e M — e m 



Q + l 



Aalog(M) (2A+(Q-l)A 2 )a 2 log 2 (M) 



O 



M 

log 3 (M) 
M 3 



2M 2 



(91b) 



Repeating the analysis of Appendix [D] with the general form 



j3(M) = X 



yields 



log(M) , ,> g 2 (M) 
M 2 



M 



Y- 



O 



log 3 (M) 
M 3 



(92) 



D[qs(t),q° s (t)} < 



Q\ 2 a 2 log 2 (M) Q 



A^log^M))-*-^ 

1-(T 

if J > Ji 



§ (Aolog(M)- 
i 



M 1_a 



o 



/log! 



' A 2 a 2 log 2 (M) \ 
4M 1 -" J + 



Q 



(Aalog(M)- 



q+i 



A 2 a 2 log 2 (M))^4 Tr 

l-C 2 

if J < J\- 



M l-A 



o 



log 3 (M) 

M 2-\a 



(93) 



Because of the M l ~ a in the denominator of the leading order 
terms, the bound for Jo > Ji converges to zero as M —> oo 
if a < 1 and diverges if a > 1. In the case J < Ji, where 
A > 1, the bound converges if a < 1/A and diverges if a > 
1/A. " " □ 



of a is larger in this case. 



B. Trade-off between e and M 

Another important trade-off is between the strength of the 
measurements e and the number of measurements required for 
effective convergence. 

Lemma 9. Allowing only e and M to vary, the distance upper 
bound (|88b|> converges to zero as M — oo provided 



£ > 



M a-l/2 



(94) 



for some a > 0. 



Proof. Consider the large M, small e regime, where the bound 
( |88b| i is approximately Bi/M. Then, since ( q /(l — C 9 ) is 
large when e is small [recall that ( = l/cosh(e)], we can 
neglect the constant first term in Eq. ((89|, 



e-^(l + rJ m )«I^ = l-l 
2Bi v ; C q C q 



3(1-0, 



so that 



f « 1 _ e r J »^L (1 + rJm)i 

and therefore, since (al - e 2 /2 in this regime, 
1 1 



(95) 



(96) 



Bi 



e Tj ™QrJ x {l + TJ m )/q 



M e 2 M 



< Ar 2a , 

where in the last line we used Eq. d9£ 



(97) 

□ 



This quantifies the previously mentioned trade-off that the 
weaker the measurements, the greater M must be to compen- 
sate, i.e., the slower the convergence rate as M — > 00. Con- 
versely, we can interpret Eq. (j94j» as saying that the measure- 
ment strength may not decline faster than the inverse square 
root of the number of measurements. 



The closer a is to the cutoff value of either 1 or 1/A, the 
slower the convergence of the bound becomes. We can con- 
clude that Eq. ( |9"0| l is a sufficient condition for convergence 
of the upper bound. Divergence of the upper bound does not 
preclude the distance D from converging to zero in the limit 
of a large number of measurements, and additional analysis 
is required to settle whether this is possible. The difference 
between the J < J\ and J > J\ cases can be understood 
intuitively in terms of the relative importance of the "pertur- 
bation" (coupling of the system to the bath, J{) and the "un- 
perturbed" evolution (Jo). It is to be expected that protection 
via the Zeno effect should be more effective when the pertur- 
bation is "weak", and as we have seen, indeed the cutoff value 



C. Fixed Nonzero Measurement Interval 

As an alternative analysis of the time-scaling issue, fix some 
At > such that the measurements are separated by an 
interval At, so that r = MAt. Under these conditions, 
what value of M minimizes the distance upper bound? Let 
6 = Jo At and A = J1/J0 and consider first the strong- 
measurement limit: 



||5||i</(M) :=e 



6M 



,Q\5 



Qe 



-xs 



M 



- 1 



(98) 
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Then, taking M to be a continuous variable for a moment, 



^ j/( M ) + l06 (£^ )WM) 

Observe that since Q = 2® - 1 > 1, 



„<5M\ 



(99) 



Qe 



-A5 



g((-l) fc + Q fc - 1 ) (rJi/M)* 



fe=2 



Q + l 



> 1, 



(100) 



the logarithm in Eq. ( |99] > is positive, as are all the other terms 
in ((99]), so df/dM > for all M > 0. Therefore on the do- 
main of positive integers M, f(M) is minimized at M = 1. 
Thus, if laboratory conditions dictate a minimum interval be- 
tween measurements, the upper bound on the distance indi- 
cates that the best strategy for minimizing the distance is not to 
let M — > oo seeking eventual convergence (which will never 
come under these assumptions), but rather to endure only one 
such interval, terminating with a single measurement event. 
The resulting bound (for M = 1) is then 



\\S\\ 1 <f(l) = e s 
QA 2 



,Q\8 



Qe 



-\s 



Q + l 

(<5 2 + <5 3 ) + 0(<5 4 ) 



- 1 



Q , 



= |(J 1 Ar) 2 (l + J Ar) + O[(J Ar) 4 ] (101) 

so that minimization of the bound can only be accomplished 
by minimizing either At or Ji . 



VI. CONCLUSION 

Two protocols have been presented for the protection from 
the environment of an arbitrary, unknown state encoded in 
some stabilizer quantum error correction (or detection) code. 
These protocols involved frequent weak non-selective mea- 
surement of either all elements of the stabilizer group of the 
code, or of just a minimal generating set. Rigorous upper 
bounds were obtained on the distance between the final state 
under these protocols and the idealized final state in the ab- 
sence of any interaction with the environment. The bounds 
demonstrate that the protocols exhibit the desired protection in 
the limit of many measurement cycles. Moreover, the bounds 
offer information about the degree of protection attainable 
with finite resources (e.g., finitely many measurement cycles), 
as well as trade-offs among the various relevant physical pa- 
rameters. 

Future research in this area may proceed along different 
lines. First, while the protocols based on non-selective mea- 
surements presented herein are well-suited to protecting en- 
sembles from the environment, they may not be ideal for pro- 
tecting individual quantum systems. Further investigation is 
needed to determine whether (and in what sense) the corre- 
sponding protocols based on weak selective measurements re- 
alize a Zeno effect resulting in protection from the environ- 
ment. A result in this direction would more directly relate the 
Zeno effect to traditional quantum error correction. Addition- 
ally, it may be interesting and fruitful to expand the class of er- 
ror correcting codes on which the protocols are based, perhaps 
to include non-abelian codes. For example, a Zeno protocol 
based on the Bacon-Shor code |41| may have advantages in 
that the measurements need be only 2-local. 



This result does not mean that the Zeno effect fails to pro- 
vide protection beyond M — 1, but rather that (our bound on) 
the protection quality gradually declines as the elapsed time 
grows. The relevant yardstick for comparison is then the dis- 
tance between the ideal and actual state in the absence of any 
protection. This distance is easily estimated using first order 
perturbation theory (the Dyson series) to be 0( J\ At). Since 
Eq. (jTOTJ shows that a single measurement already modi- 
fies the distance to 0[(JiAt) 2 ), protection is achieved pro- 
vided JiAr < 1. Subsequent measurements, or longer evo- 
lution times in the case without measurement, modify these 
estimates to 0[(J\M At) 2 ] and 0(J\M At), respectively, so 
that the conclusion about the possibility of an advantage from 
measurements with finite and fixed At are unchanged. 
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Appendix A: Isotypical Decompositions Induced by the 
Stabilizer Group 

Lemma 10. Let .M('H) denote the space of Hermitian oper- 
ators onfL = fLs ® Hb- The stabilizer group acts on A4(fL) 
by conjugation [i.e. the adjoint action Ads (A) ^ SAS], 
which describes a representation of S (~ )■ There is then 
a unique isotypical decomposition of M(fi) into subspaces 
W 9 : 

M«) = ©^=0^f s , (AD 

where each W g is an invariant subspace of the representation, 
comprising b g copies of the one-dimensional irreducible rep- 
resentation W g , such that for any Hermitian matrix A G W g , 
SAS = (-1)_ CT « {s) A. Since each W g is one-dimensional (be- 
cause S ~ I*® ™ abelian), the exponent b g is the dimension of 
the subspace W g . With Tr(l) = 2" and all other elements of 
S traceless, b g = 2 2 "- < 3(dim(H B )) 2 = 2 2k+ ® (dim(H B )) 2 
for all g. The 2 2k+< ° 1 (dim('HB)) 2 '-dimensional subspaces W g 
are all orthogonal under any Ad(S ® \ B )- invariant inner 
product on W. 

Proof. By OH Corollary 2.16], 

h * = i E Tr(Ad s )(-l)^ (5) = 2 2n -Q(dMH B )) 2 
ses 

= 2 2k+Q (dim(H B )) 2 (A2) 

under the assumption that only S = 1 has non-zero trace 
andTr(l) = dim("H) = 2 n dim(H B ), so that Tr(Adi) = 
[Tr(l)] 2 = 2 2 "dim("H B ) 2 . 

If A g £ W g and A h £ W h , then, using the Ad(S)- 
invariance of the inner product (invariance with respect to con- 
jugation by S), we get 

(A g ,A h ) = (SA g S,SA h S) = (-iy^+^{A g ,A h ) 
= (-iy^ s \A g ,A h ) (A3) 

for any S G S. If g ^ h, Lemma [TJ shows that there exists 
an S E S such that a g h(S) = 1. It must therefore hold that 
(A g , Ah) = 0, so W g and Wh are orthogonal subspaces. □ 

Appendix B: Solving Linear Recurrences 

In this appendix, we review a simple method for solving a 
class of second order inhomogeneous linear recurrences with 
constant coefficients. Let fk — di/fc-i — a 2 fk-2 = 9k be 
the linear recurrence with inhomogeneity gt — bc k for real 
numbers oi, a%, b, and c, and given initial conditions /o 
and fx. Define the characteristic polynomial to be p(x) — 
x 2 — a\x — a.2- The class of problems we will consider are 
those for which the characteristic polynomial admits two dis- 
tinct roots, i.e., where a 2 + Aci2 > 0, and where c is different 
from both roots. As in the theory of linear differential equa- 
tions, we solve such problems by finding all solutions to the 
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homogeneous problem /& — a\fk-i — 0-2 fk-2 = 0, then seek- 
ing a particular solution for the inhomogeneous problem, and 
finally identifying within the resultant affine space of solu- 
tions, the unique solution that satisfies the initial conditions. 
To that end, let 



/fe+i 
fk 



and A = 



a\ a 2 
1 



(Bl) 



Then the homogeneous linear recurrence may be written 
w k = Aw k ^ 1 = A k w°. It may be noted that the characteris- 
tic polynomial of the matrix A is identical to the characteristic 
polynomial for the recurrence, p(x), defined above. The as- 
sumption that p have distinct roots, then implies that A has 
distinct eigenvalues A±, with associated eigenvectors {P\ It 
follows that w° can be decomposed as w° = av~ 



yielding w k = A k W° 



= av + + f3v , 
= a\\v+ + p\ k _v-. Then f k = 

so the solutions to the homogeneous 

fc 



recurrence are of the form fk = 7+A+ + 7_A*i for some co 



•j k = av2 Xk 
ce 

efficients j± 

For the particular solution to the inhomogeneous problem, 
let us postulate an ansatz of fk = b'c k . Plugging this into the 
recurrence and dividing by c fc ~ 2 yields b'(c 2 — a^c — a-i) ~ 
be 2 . Since by assumption c is not a root of the character- 
istic polynomial, c 2 — aic — 0,2 7^ 0, so fk = b'c k with 
b' = bc 2 /(c 2 — aic — 02) is a particular solution to the in- 
homogeneous problem. Then the affine space of all solutions 
to the inhomogeneous problem is given by 



fk = 7 + X k + +l-X k l 



bc k + 2 



(B2) 



cr — a\c — a,2 

for coefficients j± . These two coefficients may now be deter- 



mined from the given initial conditions /o and f\ by solving 
the system of equations 



/o = 7+ + 7- + 



be 2 



— a\c — a,2 



h = 1+K +7-A- + 
which may be rewritten 



" 1 l' 




7+ 




fa C 2 


A+ A_ 




7-. 




h C 2 



be 3 



— a\c— 0,2 
be 3 



(B3a) 
(B3b) 

(B4) 



yielding 



7± = ±- 



/1 - A T /o - (c- A T ) 



6c 2 



<X\C - 



a.2, 
(B5) 



Appendix C: Inhomogeneity of the Linear Recurrence 



We seek a simple form for the inhomogeneous term in the 
linear recurrence for the sum <fi(M) ( [72) . This requires sum- 
ming both sides of |74} and using the fact that the right hand 
side leads to telescoping sums: Y^k=i^-kX(k) = X(n + 
1) — X(l). Since the limits of summation are different for 
(j)(M - 2), (f)(M - 1), and <f>{M), we will sum (|74j» over 
r/, r = 1, . . . , M and u = 1, . . . , M — 2, and add in the re- 
maining pieces of <fr(M — 1) and <fi(M) separately. 
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+ Ql3)^(M - 2) - (1 + /? + (1 + Q/3)C)0(M - 1) + 0(M) 

M-2 M M 

= EEE [£(l + + Q0)*(M-2,u,T?,r)- (1 + /3 + (1 + Q/3)^)$(M- l.u.rj.r) + *(M, u.jj.r) 



1 r;— 1 r— 1 

M 



^ [ - (1 + /3 + (1 + Q/9)£)*(M - 1, M - 1, r?, r) + $(M, M - 1, »y, r) + *(M, M, r?, r) 

ry.r— 1 
M M 

E E [/3$(M - 1, M - 1, r?, r + 1) + $(M - 1, M - 1, r? + 1, r + 1) - $(M, M - 1, t] + 1, r + 1) 

- -=i 

M M-2 

EE[( 1 + - 2 - «. V, 1) - (1 + P + OHM - 1, u, /?, 1) + $(M, «, r?, 1) 

77— 1 M— 1 

M M 

EE[~( 1 + ' 8 + ( 1 + <5^)0 $ ( M - 1, M - 1, »7, r) + *(M, M - 1, r?, r) + $(M, M, ij, r) 



(Cla) 



r/— 1 r— 1 

M M-2 



(Clb) 



77=1 r=l 
M-2 M r 

u=l r/— 1 

M-l 



+ 



EVr-W M _f)l + [- + (i + qm) E ^- 1 c M - 1 Q r? ( M _ 2 ) 



,=1 
M - 1 



+ E^- i r- 1 ^- 1 r ? (^_ 2 2 )+E^-w v „_ , 

r;=2 ^ ' ' ?7=1 ^ ' 

M-2 

= Q^2 0(1 + f3) M - u - 2 C(l + 13-0 + QK M ^ = W + f3) M - 2 t 



(Clc) 
(Cld) 



u=l 



Expand f3 = f3(M) as 



Appendix D: Asymptotic Analysis of the Distance Bound 



for some constants X > and Y > 0. Then 



1 + /? + (1 + Qflg = (1 + + (1 + QQP = (1 + + (1 + A f )X + (1 + )F + Q 

i + + Qflg = (i - + (i - QQP = a - + (1 ~m )X + (1 + o 



+ (1 + Q/3)^) 2 + 4Q£/3 2 = (1 - 



2 , 2(1-0(1-00^ , (i + Q0 2 ^ 2 + 2(i-0(i-QO^ 



M 



+ 



M 2 



+ 



Then using the fact that 



yja + x = ^fa - 



x x 
Va" ~ 8a 3 / 2 



0(x 3 ) 



(Dl) 

(D2a) 
(D2b) 
(D2c) 

J-V 

M 3 J 
(D2d) 

(D3) 
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we can compute 



l + /3-(l + Q/3)e) +4QC/3 2 =(1-0 

, X 



x(i - QQ 

M 
Y + QX 



2 Z 



Y{l-Q^)+2QX 2 
1 



(1-0J M 2+ ° \ \P ) 



(i-0 

1 ^ 



M 2 + ° ( M 3 



grc - qx 



2 £ 



(i-0 



M 2 + ° (m 3 ) 



which yields 



7+ (7+ - 7-) = 7+ ^(1 + ~P (1 + Qffl? + 4Qg/? 2 ' = (1 - + — j| + ^ + Q ^ ' 



M 2 



7- (7- - 7+) = -i-^(i+p-(i + Qm) 2 +m^= -ai - o - X ^ (q+ M 2Qg) + ° 



And using the expansion 



it may be seen that 



11 x 

a + x a a 2 



0(x 3 



Moreover, 



1 = 1 X(2-(Q + 1)0 f 1 

7+(7+-7-) l-€ " (1-0 2 M \M 2 

1 1 X(Q + 1-2QQ / 1 

7-(7-~7+) €(1-0 C(1-0 2 M + \M* 



(l + /3)[(l + /3)- 7+ ]=o(-^) 

(i^)[(i^)-7-] = (i-0 + X(2 ^™ + o( J A 2 ) 



so, assembling the pieces, we get 



7+ (1-0^ W 

§± = __0xt_ + o (J_) 
7 _ (i-0 M V M V 



(D4a) 
(D4b) 
(D4c) 

(D5a) 
(D5b) 

(D6) 

(D7a) 
(D7b) 

(D8a) 
(D8b) 
(D8c) 
(D8d) 

(D9a) 
(D9b) 



Turning now to the expressions 7 ± and (1 + /3) , observe that 



Then, 



and therefore 



Then since 



and 



f3(M) 



it follows that 



_2£ 

e m = 1 



XX 2 
M + 2M 2 + 



7 + e M = 1 + 



7 f e- X = 1 



QX -1 

7_e « £ = 1 + 



= 1 + 



(l + /3)e-w = 1 + 



(l + /3) M e- A =1 + 



y + x 2 

y + x 2 

Qy-x 2 

Qy-x 2 
x 2 



(m 3 ) 
1 

i-e 2 
Qt 1 



1 

M 2 

1 

M " 



7+ = e 



x 



(1 + /3) M 



7+ 
7- 



= 



y 



y 
y 

1- 

£M+1 

M 



2 

X 2 



1 

M 2 

1 

M 



i-e 2 

g Q 2 
1-42 

1 

1 

M 2 J 



+ 

+ 
1 

M 2 
1 

M 



m 3 y 
1 

m~ 2 



o 
-o 



1 

M 3 

m 2 y 



+ 
+ 



Y + Q(X + X 2 ) 



2^ C 



_ X 2 







1-1 


+ o( 


M J 





1 

M 



m 2 y 



m{[ QeX{x+x2) ^i 



M 



+ 



1 

M 



(m 2 ) 



m 2 y 



1 _ m y 

is ~ ~x ~ x 2 



+ 




■QJi 

e m 



e m 



) 



, 7l t 2 (2J J 1 + J 1 2 (Q-1)) 



M 



+ 



2M 2 



r g (M) 

/3(M) 



Jo > Jl 
Jo < Ji 
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Recalling that S, — C q the upper bound on the "weak" term is given by 



|W||i <r ff (iw>(M) 
= r g (M) 

/?(M) 



/3(M)A+(M) 7 f -^M) + /3(M)A_(M) 7 ^- 1 (M) - (1 + /3(M)) 



Qe rJo (rJi+rVoJx) ^ M + O (j^) J > Ji 



,1/ 



and the upper bound 



|£||i < [1 + Ti(t/M)] M -e Tj ° 



t z QJ{ 



2 1 



M 



O 



( 



1 



(D17a) 
(D17b) 



Appendix E: Correlation Functions, Spectral Densities, and 
Bath Norms 



Consider the Hamiltonians Hsb — ^2S a ® B a and Hb = 
I ® Bq. The pure-bath unitary evolution operator is Usit) = 
exp(-itBo). The bath-interaction picture bath operators are 



B a (t) = U B (t)B a (0)ul(t), B a (0) = B a 
Equivalently 

B a = -i{B ,B a (t)], B a (0) = B a 
Bt = -i[B ,Bt(t)}, Bt(0) = Bl 

It follows that after differentiating n times wrt t: 

d" 



dt n 



Bi(t) = (-ir[ n B ,BUt)], 



(El) 

(E2a) 
(E2b) 

(E3) 



where [ n B, A] := [B, [ n -iB, A]], [ B,A] := A. 
Now consider the correlation function 

(Bl(t)B ) :=Tr[e B Bl(t)B l3 ], 



(E4) 



where qb is the initial bath state. The bath spectral density is 



(D16a) 



(D16b) 



■h < Ji 



I 

the Fourier transform 
so that 



1 

2^ 



e^(Bl(t)B f3 )dt, 



t S a p{oj)6u. 



(E5) 



(E6) 



Differentiating both sides n times yields 



([ n B ,Bl(t)]B p ) = 
and in particular at t = 0: 



uj n e- iu]t 



S aP (w) dw, (E7) 



/OO 
u n S a p(u) = ([ n B , Bl)Bp). (E8) 
-co 

This last result allows us to bound the spectral densities 
in terms of the bath operator norms. One can show that 
I (AB) I < 1 1 A 1 1 ||-B|j, where the norm is the sup-operator norm 
ll42l Appendix D]. Thus 



u n S a p{ui) du> 



\([ n B ,Bi}Bp 



(E9a) 



< ||[„S 0)J Bt]||||^|| (E9b) 

< (2||S ||)"|| J Bt|||| jB/3 ||. (E9c) 

Thus, as long as ||-Bo|| and \\B a \\ are finite for all a, the 
bath spectral densities S a p(uj) must decay faster than any 
rational function as —> 00. Conversely, divergence of 
any of the moments of the spectral density implies the di- 
vergence of at least one of the bath operators. For exam- 
ple, all the n > 1 moments of a Lorentzian spectral density 
S(ui) — ( CJ _ a ,J)2 +7 2 diverge. This spectral density arises 
from exponentially decaying correlation functions, i.e., its 
Fourier transform is F^(t) oc exp(— itujQ — j\t\). 



